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RESUMEN 

Recuperar la Historia de Formacion Estelar (SFH) do una galaxia a partir de su espectro integrado es el 
objetivo central dc la smtesis de poblacion estelar. Recientes avances en modelos de sintesis evolutiva han dado 
un nuevo aliento a este viejo campo de investigacion. Modernas tecnicas de smtesis espectral incorporando 
estos avances hoy permiten el ajuste del espectro de una galaxia A-por-A. Estos ajustes detallados son utiles 
para un niimero de estudios, come h'neas de emision, cinematica estelar, y, por supuesto, evolucion de galaxias. 
Aplicaciones de esta abordaje semi-empirica a mega bases de datos nos esta enseiiando mucho acerca de la vida 
dc galaxias. El codigo de smtesis espectral STARLIGHT cs una de las herramientas la cual permite explotar esa 
combinacion favorable de datos y modelos. Como una ilustracion, mostramos como la SFHs varian para galaxias 
en distintas regiones de un diagrama de diagnostico clasico. Se detectan fuertes tendencias sistematicas a lo 
largo de las secuencias de formacion estelar y galaxias activas. Experimentos con nuevas versiones de modelos 
de sintesis son brevemente descritos. Por fin, anunciamos la disponibilidad publica tanto de STARLIGHT como 
una base de datos con ajustes espectrales detallados y productos relacionados para mas de medio millon de 
galaxias del SDSS. Esta herramienta permite exploraciones ffsicamente inspiradas del espacio de parametros, 
proporcionando nuevas maneras de navegar por el reino de galaxias. 

ABSTRACT 

Retrieving the Star Formation History (SFH) of a galaxy out of its integrated spectrum is the central goal of 
stellar population synthesis. Recent advances in evolutionary synthesis models have given new breath to this 
old field of research. Modern spectral synthesis techniques incorporating these advances now allow the fitting 
of galaxy spectra on an A-by-A basis. These detailed fits are useful for a number of studies, like emission 
line, stellar kinematics, and specially galaxy evolution. Applications of this semi- empirical approach to mega 
data sets are teaching us a lot about the lives of galaxies. The STARLIGHT spectral synthesis code is one 
of the tools which allows one to harness this favorable combination of plentifulness of data and models. To 
illustrate this, we show how SFHs vary across classical emission line diagnostic diagrams. Systematic trends 
are present along both the star-forming and active-galaxy sequences. We also briefly describe experiments with 
new versions of evolutionary synthesis models. Last but not least, we announce the public availability of both 
STARLIGHT and a database of detailed spectral fits and related products for over half a million galaxies from 
the SDSS. This facility allows more physically inspired explorations of the parameter space than is possible in 
terms of raw observed properties, offering new ways to navigate through the realm of galaxies. 
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1. INTRODUCTION 

How do galaxies assemble their stars over time, 
i.e, how do they evolve? This broad question lies at 
the heart of a large number of astrophysical frontier- 
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problems, from the internal physics of galaxies to 
cosmological issues. Theoreticians try to answer this 
question plugging in as much physics as they can 
in their models (see Abadi's talk), while many ob- 
servers tackle it using the expanding Universe as 
a time-machine, examining how galaxy properties 
change with rcdshift (e.g., Aretxaga's talk). This 
contribution deals with what can be learnt from a 
third and independent method: Uncovering the fos- 
sil record of evolution from its imprints on galaxy 
spectra. This semi- empirical approach has become 
highly attractive in the past few years, given the 
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avalanche of data from cosmologically shallow sur- 
veys, and the enormous progress in our ability to fit 
galaxy spectra on a pixel-by-pixel basis using state- 
of-the-art evolutionary synthesis models as those re- 
viewed in Bruzual's and Coelho's contributions. 

Taking full advantage of this favorable combi- 
nation of abundance of data and models requires 
spectral synthesis tools to extract information about 
age (t), metallicity (Z) and the detailed star for- 
mation history (SFH) encoded in observed spectra. 
There are now several such tools around, differing 
in both technical and astrophysical aspects (eg, Ma- 
teu's talk). Some account for extinction and/or kine- 
matics, others don't. Some impose simple Z-t rela- 
tions or use a fixed Z, while others treat Z and t 
independently. Some model indices, others the full 
spectrum. Some fit the SFH in a non-parametric 
fashion, while others compare the data to a library 
of precomputed models. Some prefer to compress the 
input data, others the output parameters. The list 
goes on and on. . . We skip the impossible task of cov- 
ering all this ground by referring the reader to two 
reviews by our group (astroph/071899 and 071902), 
dedicated to basic aspects of spectral synthesis and 
the recent literature in the field. 

These few pages focus on results obtained with 
our STARLIGHT synthesis code. Rather than self- 
advertisement, this is done with the specific purpose 
of illustrating the sort of science doable with this 
code, specially when applied to mega data sets. The 
motivation for this is that both STARLIGHT and 
products of its application to ~ all SDSS galaxies are 
now available for public use. In www. starlight.ufsc.br 
the reader will find the code itself, A-by-A fits for 
573141 galaxies from the SDSS DR5 and a long 
and diversified list of derived properties (from stel- 
lar masses to emission line fluxes). This database is 
about to ~ double with 694135 DR6 spectra fitted 
with new evolutionary synthesis ingredients. 

We start with a quick introduction to 
STARLIGHT, its deliverables and our VO-like 
database (SJ2]), including a first ever comparison 
of results obtained with DR5 data modeled with 
the standard Bruzual & Chariot (2003) models 
with DR6 data fitted with newer models. As an 
invitation to our database, we illustrate the power of 
our detailed spectral synthesis approach by showing 
how SFHs vary as a function of location on classical 
emission line diagnostic diagrams (ij3)). 

2. THE STARLIGHT/SEAGAL PROJECT 

2.1. The code 
STARLIGHT (Cid Fernandes et al. 2005) com- 
bines Ni, spectra from a user-defined base of indi- 



vidual populations in search of linear combinations 
which match an input observed spectrum. The fitted 
coefficients define an A'^.-dimensional population vec- 
tor (light fractions at a reference A) . For SFH studies 
it is useful to use a base of instantaneous bursts of 
different t's and Z^s, but anything else can be used. 
Kinematics is allowed for, as is reddening (according 
to any law) . Papers by our Semi Empirical Analysis 
of Galaxies (SEAGal) collaboration discuss the code 
and the results of its apphcation to ^ 0.5M galaxies 
from the SDSS. For an in depth description of the 
code, its pros and cons and possible uses a detailed 
45-pages long manual is available. 
STARLIGHT itself outputs: 

• The full synthetic spectrum Mx. 

• The light-fraction population vector x. 

• A mass-fraction population vector fl (only 
meaningful for properly defined bases). 

• Stellar velocity dispersion (cr^) and shift (w+). 

• Stellar extinction {Ay). 

This is already enough for those interested in, say, 
stellar kinematics, or if all you need is a decent stellar 
template to subtract from your data to aid emission 
line work. In neither case you would care about x 
nor fl. For galaxy evolution work, on the other hand, 
it is precisely the population vector which matters, 
since this is where the SFH information is. It can 
be handled in numerous ways to produce things like 
time dependent star formation rates, SFR(t). 

2.2. The database: www.starlight.ufsc.br 

In addition to the data listed above, our 
STARLIGHT-SDSS database contains a series of 
other products: 

• Emission line properties (fluxes, equivalent 
widths, widths, S/N) for many transitions, measured 
from the starlight-subtracted spectrum. 

• Stellar masses (M^,), and mass-to-light ratios. 

• Summaries of the population vector, like mean 
stellar ages and metallicities. 

• Basic SDSS data, like coordinates, redshifts, 
magnitudes, sizes and images. 

Fig. [T] shows a screen-shot of one object of a list 
produced with a SQL selection. The page contains 
the spectrum, the STARLIGHT fit and aU that is 
known about the galaxy, as well as its SDSS thumb- 
nail! A particularly attractive feature is that one 
may search and select by any combination of physi- 
cal and observed properties. This allows one to play 
all sorts of games. The bottom panel in Fig. [T] shows 
an example. We have selected galaxies within 0.05 

^The whole interface was developed and is maintained by 
an undergrad student (W. Schoenell)! 
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Fig. 1. Top: Screen-shot showing an example database 
query. Bottom: Galaxies of same stellar mass 
{m^^-°°^°-°^ Mq) and redshift (0.07 ±0.01), low inclina- 
tion [b/a > 0.9) sorted according to the STARLIGHT- 
derived (light-weighted) mean stellar age: logt[yr] =8.3 
to 10.3, increasing to the right and from the top. 

dex of = 10"Mo, z = 0.07 ± 0.01 and produced 
a mosaic of 25 images sorted according to the mean 
stellar age. The progression from late to early type 
morphologies is evident to the eye. 

2.3. Elementary (but useful) warnings 

It is never too much to emphasize the following 
obvious, but often overlooked, points: (1) The inver- 
sion from an observed spectrum to a set of param- 
eters is far more complex and degenerate than one 
would hope. Yet, as long as one does not fall into 
the temptation of embracing over-detailed descrip- 
tions of SFHs, all modern synthesis methods seem to 



converge to the same overall result. Furthermore, av- 
eraging results over hundreds of galaxies greatly alle- 
viates uncertainties in individual SFHs (Panter et al. 
2007). (2) To go from an observed galaxy spectrum 
to a population vector and its associated SFH one 
must chose one among many possible bases. Evolu- 
tionary synthesis models themselves are undergoing 
constant revisions, approximately at the same rate 
of that of improvements in their two main ingredi- 
ents: stellar tracks and spectral libraries — this very 
same volume contains 2 such updates! (3) Finally, 
observational data are never problem-free. We have 
chosen to fit anything spectroscopically classified as 
a galaxy in the SDSS! This includes quite a bit of 
garbage, which must be filtered out with appropri- 
ate quality flags. 

2.4. New bases + new data = more work 

There are news about points 2 and 3 in the above 
list. First, we point out that all published results by 
the SEAGal collaboration rely on the 2003 version 
of the Bruzual & Chariot models with the STELIB 
library and Padova 1994 tracks applied to over 0.5M 
SDSS DR5 spectra. But there is a lot more coming. 

As discussed by G. Bruzual elsewhere in this vol- 
ume, there are newer versions of these models, differ- 
ing both in the treatment of late evolutionary phases 
and in the spectral libraries. We have carried out a 
series of experiments with some of these new models, 
and found that: (a) Not surprisingly, the new evo- 
lutionary tracks have negligible impact on our opti- 
cal spectral synthesis (they only affect the IR). (b) 
Fits with the MILES library (Sanchez-Blazquez et 
al. 2006) provide substantial improvement with re- 
spect to those obtained with STELIB, particularly 
for old systems. Besides smaller spectral residu- 
als, these new models correct some pathologies in 
the derived physical properties of passive/elliptical 
galaxies, like weird-looking distributions of sources 
in the t-Z plane and (slightly) negative values of Ay . 
Things are thus improving on the modelling front. 

There are also news about the data: The SDSS 
reduction pipeline changed, and DR6 spectra are dif- 
ferent in shape and amplitude from the DR5 ones. 

We have just finished fitting 694135 DR6 spectra 
with the new MILES-based models, and will soon in- 
corporate these data to our database. Fig. [2] presents 
a preliminary comparison of old and new fits to three 
passive and two SF galaxies. Some of the changes 
mentioned above are visible in these examples. The 
similarities and differences revealed by practical tests 
like this provide useful feedback to model makers, 
and portray a summary of the degree of maturity 
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DR5 data + BC2003 models (w/STELIB) 



DR6 data + New EC models (w/MlLES) 
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Fig. 2. Fits of the same five SDSS galaxies, but using DR5 + Bruzual & Chariot (2003) models with the STELIB library 
(left) or the DR6 version with a newer version of the models with MILES (right). Black and red are the observed and 
synthetic spectra, with masked regions plotted in magenta. The histograms show the % light fractions at 4020 A from the 
STARLIGHT decomposition onto a common set of 25 ages x 6 metallicities. Solid lines are 0.5 dex smoothed versions 
of our population vector. (Given the 4 dex dynamic range in t^, smoothing by 0.5 dex is qualitatively equivalent to 
reducing the base to 8 ages, perhaps still a bit over-optimistic.) Despite the broad-brush agreement, there are substantial 
differences both in the input data and in the fitted SFHs. Notice also the different spectral ranges. 



in the field. While these new models represent a 
perceptible improvement, before jumping to conclu- 
sions one must realize that there are many others 
available, some of which address known deficiencies 
in standard models (like a/Fc). Examining all pos- 



sibilities and digesting the results will take time and 
patience. It will be particularly interesting to inves- 
tigate the effects of different stellar evolution tracks, 
even if that requires extending the spectral domain 
beyond the optical (e.g., Eminian et al. 2008). 
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Fig. 3. BPT/ "seagull" diagram for 131421 SDSS galaxies with superimposed SFHs as derived by STARLIGHT. The 
dashed lines show the SEAGal (Stasiriska et al. 2006, in blue) and Kewley et al. (2006, orange) AGN/SF dividing lines. 
Each box shows the t-hy-t median (thick black line), 16 and 84 percentiles (thinner brown lines) of the time dependent 
SFR per unit mass for galaxies within the box. The inset shows the scale (units are yr^^ for SSFR and yr for lookback 
time t). Note the systematic progression of the ratio of past to present star formation as one moves down the SF (left) 
wing, and from the Seyfert to the LINER branches (above and below the seagull's right wing, respectively). 



3. STAR FORMATION HISTORIES ACROSS 
THE BPT DIAGRAM 
Emission lines reflect the physical conditions of 
the ISM and the nature of the ionizing source. Since 
Baldwin, Phillips & Terlevich (1981), diagrams in- 
volving pairs of line ratios are used to diagnose 



galaxies as star-forming (SF) or active (Seyfert or 
LINER), and to derive properties like the nebular 
metallicity {Z„eb), N/0 ratio and ionization param- 
eter. An interesting question to ask is: How do the 
SFHs of galaxies change across such diagnostic di- 
agrams? Given that the ISM conditions arc more 
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a consequence than a cause of galaxy evolution, this 
question is kind of an exercise in reverse engineering. 

Fig. [3] shows the classical [O III]A5007/H/3 vs. 
[NII]A6584/Hq' BPT diagram. The background in 
this plot consists of emission line measurements 
for 131421 galaxies from our STARLIGHT-SDSS 
database). The SF and AGN sequences are clearly 
visible, branching out like the wings of a seagull. We 
have chopped this seagull into small boxes and com- 
puted the statistics of the t-by-t SFH within each box 
(like in Asari ct al. 2007, but including AGN). The 
thick-black line shows the time evolution of the me- 
dian Specific SFR (SSFR), whereas the thiner-brown 
lines are the 16 and 84 percentiles. 

Systematic variations in the evolutionary pattern 
are seen all across the diagram. Along the left wing, 
for instance, all galaxies show current SF (as ex- 
pected), but the balance between current and past 
SFR changes by over an order of magnitude. That is 
why galaxies at the top left, where \ow-Zneb galaxies 
live, look much younger than the metal-rich systems 
at the bottom of the sequence. Mass increases along 
this sequence (the M-^-Zneb relation), so this pattern 
is ultimately a consequence of galaxy downsizing. 

Strong trends are also evident in the AGN wing. 
The stellar population mixture becomes increasingly 
skewed toward older ages as one moves up the AGN 
wing, in good part due to aperture effects. Even con- 
centrating only on the region far away from the body 
of the seagull (eg, beyond the Kewley line, where 
contamination by off nuclear SF regions is minimal) , 
one sees trends: Seyfert 2s (upward in the plot) are 
systematically younger than LINERs. In most AGN 
the ongoing SF is but a fraction of what it has once 
been, and this fraction tends to zero as one walks the 
BPT diagram from the Seyfert to the LINER terri- 
tories. Several other things change on the way too: 
Ay and line luminosities go down, while Mi,, Z^, and 
a/Fe increase. 

To have a closer look at this trend. Fig. [4] shows 
SFHs in the L[oiii]M^» vs. Mi, plane, where M, is 
the indirectly computed black hole mass (from cr^,). 
The y-axis discriminates between Seyferts and LIN- 
ERs (Kewley et al. 2006), and is usually interpreted 
as a proxy for the Eddington ratio (Heckman et al. 
2004). Only 16411 galaxies "pure AGN" are shown 
in this plot. Despite the relatively compressed M^, 
scale (AGN live in massive galaxies), it is interest- 
ing to note that for a fixed L[oiii]/M, the SFHs are 
essentially independent of galaxy mass. On the con- 
trary, for any given Mj, the "Specific Black Hole Ac- 
cretion Rate" (i.e., M,/M,) is obviously linked to the 
recent SFH: Black holes are growing more efficiently 
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Fig. 4. As Fig.O but in the "Eddington ratio" vs galaxy 
mass plane, and restricted to the 16141 galaxies above 
Kewley's maximal-starburst line. The inset shows the 
relation between the specific SFR over the last 10* yr 
and a proxy for the "specific black hole accretion rate" . 

in galaxies which are forming stars more efficiently. 
This is shown more directly in the inset, which plots 
the relation between L[oiii]/-A/, and the mean SSFR 
over the last 100 Myr. 

A few years ago this would be called "the 
Starburst-AGN connection". To find out whether 
this fits with negative AGN-feedback scenarios fash- 
ionable nowadays, read our next papers. Or else, do 
it yourself at www.starlight.ufsc.br! 
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